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Here we study an anisotropic model of the universe with constant energy per particle. A decaying 
cosmological constant and particle production in an adiabatic process are considered as the sources 
for the entropy. The statefinder parameters {r, s} are defined and their behaviour are analyzed 
graphically in some cases. 

PACS numbers: 

The inclusion of cosmological constant into Einstein field equation was done primarily by 
Einstein himself with a view to model a static universe. Later he discarded the cosmological 

constant (and commented [1] " the biggest blunder of my life"), on the basis of 

the observational data supporting the expansion of the universe. Subsequently, it was 
reintroduced and then omitted to match with the observational evidences. 

Recently, cosmological models with a cosmological constant are strong candidates to 
describe the dynamics of the universe. The observational results [2, 3] from type la 
Supernova suggests that our universe is accelerating. The increase of the expansion velocity 
can be accounted for by introducing a repulsive force which results a negative pressure 
in the energy-momentum tensor. Cosmological constant (or term, time dependent) is a 
good candidate to generate negative pressure and to account the vacuum contribution to 
the energy-momentum tensor. Further, the present small value of A [4, 5] suggests that 
it is not a constant but have a dynamical evolution. As universe expands, the effective 
cosmological term evolves and decreases to its present value. This is also supported 
strongly from observational point of view, particularly, the anisotropy of the cosmic mi- 
crowave background radiation, the supernova experiments and estimated age of the universe. 

Moreover, the A-term can be considered as a type of dark non-baryonic matter (or dark 
matter) in the sense that it is not gravitationally clustered at all scales as the usual matter. 
Lastly, it is to be noted that in the literature there are cosmological models with particle 
production which arc compatible with experimental evidences. 

In this work, using the adiabatic condition, the dissipative pressure is considered due to 
decaying of a cosmological term and particle production. The thermodynamical condition is 
necessary for the conformity of the cosmic microwave background radiation due to creation 
of photons. 

In thermodynamics, the energy conservation equation is 



nTda = dp — [fi + Ta)dn 



(1) 



where the chemical potential fj, satisfies the Euler equation 



iP + p) 



Ta 



(2) 



n 
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So eliminating /x between (1) and (2) we have 



r,Ta = p-^-P^n (3) 
n 

(' . ' stands for time derivative) 
The number density 'n' satisfies the continuity equation 

n + ne = * (4) 

with 6, the expansion scalar and \p, the source of particles. Now the second law of 
thermodynamics can be written as 

5^ = riCT + C7* (5) 
with 5" = ncru" as the entropy four vector. 
Now, considering the energy-momentum tensor as 

T^"" = {p + p^ )u^'vy - p^g^"' + hg^"" (6) 
the energy conservation equation (T'f^ = 0) results 

p + (p + p^)e = -^ (7) 

Here p is the matter density, is the sum of the usual thermodynamical pressure p and 
a dissipative pressure 11 (i.e., Pj, = p + H), u^^ is the four velocity vector and A = A{t), a 
varying cosmological constant. 

Now substituting a and n from equation (3) and (4) and using the energy conservation 
equation (7) we have from equation (5), the expression for the entropy as 



+ C7* (8) 



SwG n 

Considering the particle production process as an adiabatic process (i.e., a = 0) one writes 

A + to)* + ne = o (9) 

The line element for the Kantowski-Sachs (KS) model is given by 

The Einstein field equations for the matter content of the form (6) is given by 

2^ + 5 + ^ = -8'^^(^' + n) + A (10) 



^ + ^ + ^^=-87rG(p + n)+A (11) 
a b ab 
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and 

Assuming barotropic equation of state p = ep (0 < e < 1), one can eliminate p, p and 11 
using the field equations (10)-(12) and equation (9) to obtain 



with e = ^ + 2| = 3F. 

a ' b 

In the present model, particle production is possible due to variation of A and due to 
dissipative pressure. We now assume, the contribution due to these terms are proportional 
namely, 

A * 

= w{p + p)- (14) 



SttG " n 

with uj as the constant of proportionality. It is to be noted that vanishing of \1/ implies 
both A and entropy to be constant. 

Combining equations (9) and (14), the dissipative pressure is given by 

n = -(l+u;)(p+p)^ (15) 
Also taking into account relations (7), (9) and the field equations, one gets 

Now choosing A = aH^ as proposed by Viswakarma [6] in recent past, we see that the 
particle source corresponding to the A term is 



* = PnQ (17) 
where /3 is a phenomenological constant [7] . 

Using these relations, (namely, (15)-(17)) in the field equations to eliminate p, p and 11, 
we have the differential equation in the scale factors as 

4 + + - (1 + '^)/^} - - (1 + - (1 + '^)/5>] + + - (1 + '^)/^> 

~f (l + ^){l-(l+'^)/5} = (18) 
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Also the differential equation for A (i.e., eq.(14)) using equations (15) and (17) can be 
solved as 

A = .m + e){ab^) ^j{ab-) ^ ^- + 2- j (^^ + 2-- + ^ j .^i + Ao 

Now for solution let us assume a power-law form for b as 

b = bnf^ , bo, v are constants. (19) 
Then from the field equations (10) and (11) the differential equation for 'a' becomes 

t^a - v{2u - l)a + vta - ^i^^^'^^a = (20) 
For fc ^ 0, it has simple solution for i/ = 1 and \ as follows: 

/i+4 

v=l : a = ooiV , 6 = 6ot (21) 



1/=^ : a = aosinh ^^^^i^ j , b = bot^ (22) 
However, for general the solution can be written as 

a = I ao 7(1-3.) [TT^^i'-l + ai 7 d-a., [jT^^t'-'] \ (23) 
2TT^ {l-i^)bo -2U^'{1- iy)bo 'j 

where / is the usual modified Bessel function. 

Note that for the solution with i/ = 1 the constraint satisfied by the constants appearing 
in equation (18) has the simple form 

(l + e){l-(l + u;)/3}(3-a+|) =2 (24) 
Further, for k = 0, equation (20) has the general solution 

a = aof-" + ait^-^" (25) 

In recent past, Sahni et al [8] proposed statefinder parameters to discriminate between 
different dark energy models. In FRW model with scale factor a{t), the two parameters are 
defined as 

d r — 1 

with H and q as the Hubble parameter and deceleration parameter respectively. These 
dimensionless parameters can be written in general form for any space-time as 

3H H r~l H 
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In the recent model we shall determine these statefinder parameters for different solutions 
and examine their behaviour. 

For 1^ = 1, k ^ 0, q, r and s are all constants and A a t~'^ with expressions: 

1-e , 9g 6 I k Ao 

^=^'^ = '-(^'^ = ^'^ = V 

Note that s is positive definite for all k while the universe will be accelerating for closed 

model and decelerating for open model. For other or for flat model the expressions of 
these parameters are very complicated so we have shown their behaviour graphically only 
for 1/ = i with fc ^ or A; = 0. 

In figures 1 and 2, we have shown the time variation of q and r-s curve respectively for 
ly = ^, k 0. The figure for q shows that the universe starts from an decelerating phase 
(radiation era) to an accelerating phase and finally becomes a ACDM model. This is also 
reflected in the r-s curve. The branch on the r.h.s. of the asymptote (in Fig. 2) represents 
the universe from radiation era to dust phase (when s becomes infinity) while the left 
portion of the curve continues from the dust phase and goes gradually to ACDM model. 

For flat case (i.e., fc = 0) with = ^ the time variation of q and r-s curve arc shown 
in figures 3 and 4. These are not of much interest from observational point of view as the 
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model does not represent any accelerating phase of the universe (since q is never negative) . 

Therefore, the present model with non-flat universe and v = ^ can describe the evolution 
of the universe from radiation era to ACDM model but beyond radiation in the past is not 
describable by the model — probably some quantum theory is necessary to describe the 
evolution of the early era. 
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